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Linear induction motors (LIMs) make performing a direct linear motion 
possible without any mechanical rotary to linear motion transforming parts. 
Obtaining a precise mathematical model of such type of motors presents a 
difficulty due to time varying parameters and external load disturbance. This 
paper proposes an adaptive backstepping controller structure based on 
lyapunov stability for controlling a LIM position. Which can guarantee the 
annulment of position tracking error, despite of parameter uncertainties. 
Parameter update laws are extracted to estimate mover mass, friction 
coefficient and load force disturbance, which are assumed to be constant 
parameters; as a result, compensating their undesirable effect on control 
design. Then, load disturbance estimate is replaced with an artificial neural 
network (ANN) to reduce the estimation error. The numerical validation has 
shown better performance compared to the conventional backstepping 
controller, and proved the robustness of the proposed adaptive controller 


design against parameter changes. 
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1. INTRODUCTION 

Linear induction motors (LIMs) are widely used in industry and transportation applications [1], [2]. 
This type of motors offers the ability to perform a direct linear motion without the need of any mechanical 
rotary to linear motion transforming parts [3]. The driving principle of the LIM is similar to that of a rotary 
induction motor (RIM), however, the control characteristics are more complicated: operating conditions such 
as speed of mover, temperature and rail configuration result variations in LIM parameters [4], [5]. Field 
oriented control (FOC) is considered as the most used high performance control method for AC motors. Its 
main idea is to make the control design similar to that in a separately excited DC machine, by decoupling the 
dynamics of thrust force and rotor flux [6]-[8]. Unfortunately, a good performance of such control can only 
be achieved by using a precise plant model, which is very difficult, because unlike RIM the model of LIM 
presents more complexity due to parameter variations. Many researches in modeling LIM dynamics have 
been done, taking into consideration parameter variations and end effect phenomenon [3], [6], [9]-[11]. 
However, modeling uncertainties are always present due to undesirable behaviors and load disturbance. 

The advancement in nonlinear control field has offered multiple solutions to deal with modeling 
changes. One of these solutions is backstepping control. This nonlinear control technique has been widely 
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applied in controlling LIMs [12]-[16]. The design idea is the decomposition of a complex high order control 
structure into low order steps. Where in every step a new reference is chosen to stabilize the previous one 
using “virtual control” variables, until reaching a real control input that garentees the stability of the global 
structure using lyapunov function. Moreover, using the adaptive version of backstepping control insures 
attaining the global structure’s stability despite of parameter variations, thanks to parameter update laws [12]- 
[16]. The obtained results in previous works have shown better performance comparing to FOC when facing 
modeling uncertainties [16]. However, the linearity assumptions in the design may not always be true. The 
load disturbance for example is neither linear nor constant in most cases. 

Parallel to fast development in robust and adaptive control techniques, artificial neural networks 
(ANNs) have been widely applied to system identification, and identification based control. They are known 
to be universal approximators. ANN based controllers have been proposed for a variety of control 
applications that can provide closed-loop stability [17]-[21]. In this paper, an adaptive backstepping control 
structure is proposed in order to track a desired position trajectory of a LIM, under the assumption of 
unknown moving part mass, friction coefficient and load disturbance. Then, an ANN estimation structure is 
introduced to replace the external load disturbance estimate. The paper is organized as follows: section 2 
introduces the proposed adaptive backstepping controller design with ANN estimator. Section 3 presents the 
numerical validation results followed by discussion. Finally, a conclusion is summarized in section 4. 


2. ADAPTIVE BACKSTEPPING CONTROL 

The construction principle of a LIM is shown in Figure 1. The primary can be seen as a “cut-open 
and rolled-flat” RIM primary, and the secondary usually consists of a sheet conductor, with an iron back for 
the return path of the magnetic flux. For a neglected relative velocity, the LIM primary can be considered as 
infinite, as a result, the end effects may be ignored [13]. 


Entry rail Secondary sheet Secondary back Exit rail 
eddy current iron eddy current 


Figure 1. LIM construction 


The mathematical model of the LIM in the d-q axis reference frame can be modified from a three 
phase Y-connected RIM, which gives [12], [13]: 


tte (B+ (2) R) ia teta t tart i yb seh 0 
He = Estas (Ro + (BY fe) bp Zt ben 0 tbe + EV O 
Her tata paid, 6) 
Star w= tate iag -ECV — Pav). bar — Oar e 
F, = Ky(barigs — Pqrias) = MZ + fov +F, (5) 


where R, is the stator resistance, R, is the rotor resistance, L, is the stator inductance, L, is the rotor 


inductance, L is the mutual inductance, P is the pole pairs number, h is the pole pitch, ø = 1 — (L?,/L,L,) 
is the leakage coefficient, V, is the linear synchronous velocity, v is the linear mover velocity, ig, is the 
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d-axis stator current, iqs is the q-axis stator current, ar the d-axis rotor flux, @qr is the q-axis rotor flux, 


Vasis the d-axis stator voltage, V,, is the q-axis stator voltage, 


Ky = 3PrLm/(2L,h) is the force constant, Fo 


is the applied electromagnetic force, F, is the load force disturbance, M is the total mover mass and fis the 
friction coefficient. The basic backstepping design idea is the decomposition of a high order nonlinear system 
into small first order subsystems, using “virtual control” variables [22]. Where every control variable is 
chosen in a way that can stabilize the previous subsystem, introducing a corresponding error that can be 
compensated by selecting a convenient control input via lyapunov stability [7], [13], [14], [23]. 


2.1. Controller design 
Introducing the position tracking error: 


e1 = dres —d 

its dynamic gives. 
ey ee ar eee 

We introduce the first lyapunov function candidate, 
vV = zef 

it’s time derivate gives, 
A = e1] = e1(dref —v) 

if the first stabilising function is chosen as, 
Vref = kye, + dref 

where k is a positive constant, the dynamic of V; gives, 
V, = —k,e? < 0 


which insures the stabilisation of the closed-loop control. 


(10) 


(1) 


(12) 


(13) 


(14) 


(15) 


However, v is a “virtual control” input, a new tracking error is defined in (16). 


e2 = Vref —V 


(16) 


To ensure the elimination of the position tracking error despite of modeling inaccuracy and load 
disturbance, an integral action is introduced to the velocity reference as follows [15]: 


Vref = kye, + dref + k,'e,' 


where €,' = f e,(t)dt and k,' is a positive design constant. 
Deriving the velocity tracking error to time gives, 


€z = Vres — Ù 
= dre, + kı (die, - v) +k'ie, -Ż+Év +2 
so the dynamic of e, can be rewritten as, 
é = —kye, + ky'e,' + ez 
define a new lyapunov function, 


1 1 ky! 
V, = zei +>e are 
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deriving to time and using (18) and (19) results, 
V, = —k,e? aa k2e2 + e2 [e (1 a ke + kı”) + ki kiej“ +e(k, + k2) + dref ga Z4 D. v + L] (21) 


where k, is a positive constant, D= f,/M is the normalized friction coefficient and L= F, /M is the 
normalized Load disturbance. Now, the force control input F,_+ef that guarantees achieving the desired 


position trajectory is to be chosen. If we choose (22), 

Fe ref = M[e (1 — k? + k1") + ky'kyey! + ez(k1 + k2) + d,ef +D. v + L] (22) 
we get, 

V, = =ke? — kze? < 0 (23) 


which ensures the stabilization of the closed-loop control. 


A 


However, M , D and L are considered as unknown parameters, so their estimates M , D and L 
are used instead, (23) becomes (24). 


Fe ref = Mle (1 — k? + ky!) + ki'kie,' + ez(k1 + k2) + dre +D. v + L] (24) 


Next step is to extract estimation update laws. Define the following parameter estimation errors. 


M=M-M (25) 
D=D-D (26) 
L=L-1L (27) 


Using these definitions and the applied force control (24), the dynamic of e, in (18) gives, 

; = > im 

êz = —€; — kz€62 + D.v +L +3Ê (28) 
where B = [e,(1 — k? + ky) + ky'kyey' + e2(ky + k2) + dref +D. v + Î]. 


Finally, the system’s final lyapunov function that contains all error signals is build, in order to 
extract update laws of the estimated parameters, 


—1,2,1,2,%', 2 1 7241724172 
V3 er +5 en + al tzau +z P +558 > 


where ôi, ô and 63 are positive constants. 


Assuming M , D and L are constant parameters, It can be written, 


M = -M (30) 
Ď = -Î (31) 
Leer (32) 


Using (30), (31), (32) and the error dynamics in (19) and (28) we calculate the time derivative of 
(29). After substitution and simplification, the result is: 


V; = —k,e? kzež + "(eB =h) + D (ev - ZD) +1 (e, - Z2) (33) 
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The chosen update laws to render the system’s global lyapunov function negative are as follows: 


M = 6,e.8 (34) 
D = 6,e,v (35) 
L = d3e, (36) 


The resulting lyapunov function dynamic gives, 
Vz = —k,e? — ke2 < 0 (37) 


which guarantees stabilising of the global system and achieving desired position trajectory, despite the 
presence of parameter uncertainties and load disturbance. 


2.2. ANN load disturbance estimator 
From (36), the load disturbance estimate can be written as (38). 


L=63fe dt (38) 


In the next approach, we will replace the load disturbance estimate with an ANN structure in order 
to reduce the load disturbance estimation error. We have chosen an ANN structure with single input e, and 
single outputLyy. The chosen ANN contains three layers: the input layer has one neuron and no activation 
function, the hidden layer consists of 5 neurons with tangent-sigmoid activation function, and the output 
layer has one neuron with purelin activation function. The selected ANN was trained off-line using previous 
controller results. 

The basic feed-forward three-layer neural network model is given as [24], [25]. 


Ye (XW) = fE Wie j 9 (Lies Wii Xi + bi) + bj) (39) 


where yjy,..... ,Y~ are output variables, f is the output layer activation function, x ,....,X, are input 
variables, wę j; and b j are respectively hidden to output layer weights and bias, g is the hidden layer 


activation function, w; j and b; are respectively input to hidden layer weights and bias. 


Applying the the neural network form in (39) to the load disturbance estimate in (38), our ANN 
structure can be written as (40). 


Evy = fÈ- W2i9 (w, je. + bi) + b;) (40) 


The final adaptive backstepping controller design with ANN estimate is summarized in diagram of Figure 2. 


Adaptive 
backstepping 


Figure 2. Adaptive backstepping control with ANN estimator 
3. RESULTS AND DISCUSSION 
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In the following, we run numerical simulation in which the performance of the adaptive 
backstepping controller under different operating conditions is tested (Mass, friction and load disturbance 
variations). For that purpose, we run numerical simulation for the following cases: 

— Case (1): variation of mass equals to 3xM. The results are shown in Figure 3. 

— Case (2): variation of friction coefficient equals to 10xD. The results are shown in Figure 4. 

— Case (3): a force disturbance equals to 20N is injected between t=5s and t=7s. The results are shown in 
Figure 5. 

The obtained results are compared to conventional backstepping controller performance. The 
machine parameters used for numerical validation are presented in Table 1. 


Table 1. LIM parameters 


Parameter Value 

Rs 3.4 (Q) 
R; 1.95 (Q) 
Ls 0.1078 (H) 
Lr 0.1078 (H) 
Li 0.1042 (H) 
P 2 

MM 5.47 (Kg) 
fe 2.36 (Nm.s/rd) 


The chosen trajectory is a periodic step of +0.1m and the controller’s parameters and update laws are selected 
as: k4 = 10, ki = 0.02, kz = 120, 6, = 0.001, 6, = 0.01, 63 = 500. 


3.1. Mover mass variation 

In the first case, the results show that the adaptive backstepping control is insensitive against mover 
mass change. Even when the mass value rose to triple, the controller is always insensitive and can track the 
desired position trajectory successfuly with fast response time and no overshoot Figure 3a. In the other hand, 
the conventional controller has suffered a tracking response overshoot, before it could achieve the reference. 
As a result, there was a delay in response time. The electromagnetic force, stator current and rotor fluxes are 
shown in Figures 3(b)-(d) consecutively. 


“Position reference 600 
— Adaptive backstepping 
-- Backstepping 400 

7 


200 


0 


Fe (N) 


-200 


-400 


-600 


Time (s) 


Time (s) 


(a) (b) 


Tgs (A) 


Time (s) Tir 


(c) (d) 


Figure 3. Adaptive backstepping control of LIM position with 3xM, (a) position tracking response, 
(b) applied electromagnetic force, (c) stator current Iqs, and (d) rotor fluxes 


3.2. Friction coefficient variation 

In case 2, where the friction coefficient value has risen by 10 times, both controllers had nearly no 
noticeable effect. The response time is fast and there is no overshooting as shown in Figure 4(a). That shows 
that backsteppnig control is insensitive against friction variations too. The electromagnetic force, stator 
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current and rotor fluxes are shown in Figures 4(b)-(d) consecutively 


Position reference 
0.15 —— Adaptive backstepping H 
-- Backstepping | 


Fe (N) 


Time (s) 


(b) 


(c) (d) 
Figure 4. Adaptive backstepping control of LIM position with 10xD: (a) position tracking response, (b) 
applied electromagnetic force, (c) stator current Iqs, and (d) rotor fluxes 


3.3. Load disturbance 

In case 3, viewing Figure 5 we can see the effect of load disturbance on the conventional controller, a 
non-null static error was generated causing a failure in following the position trajectory as seen in Figure 5(a). 
In the other hand, the adaptive backstepping controller was able to compensate the effect of load disturbance 
in fast time. In addition, using the ANN estimate has made the controller insensitive against load disturbance, 
making the position tracking successful with no noticeable effect. The applied electromagnetic force, stator 
current and rotor fluxes are shown in Figures 5(b)-(d) consecutively. 


Position reference 


| 0.1 ee 
0.154 Backstepping 19.09. . Sate | 600 
— Adaptive with ANN 10.08 r 
---- Adaptive without ANN 1.48 5.8 400 


d(m) 


Fe (N) 


0 1 2 3 4 5 6 7 8 0 1 2 3 4 5 6 7 8 
Time (s) Time (s) 


(a) (b) 


lgs (A) 


Qdr , Qqr (Wb) 


Time (s) Time (s) 


(c) (d) 


Figure 5. Adaptive backstepping control of LIM position with applied load disturbance equals to 20N; 
(a) position tracking response, (b) applied electromagnetic force, (c) stator current Iqs, and (d) rotor fluxes 
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The previous simulation results show great performance and prove the robustness of the proposed 
adaptive backstepping control structure, thanks to parameter update lows using lyapunov function. Which 
insures the global system’s stability despite of parameter uncertainties. The parameter estimates are presented 
in Figure 6. We can observe that the estimation errors are null for mover mass in Figure 6(a) and friction 
coefficient in Figure 6(b). For the load disturbance estimates in Figure 6(c), we can see that the estimation 
error was considerable and can’t be neglected for the structure with no use of ANN. This problem is the main 
reason we proposed the replacement of the load disturbance estimate with an ANN estimator. This structure 
has reduced the estimation error. It can be noticed that load disturbance estimation converges to the real 
value, with better and more acceptable response time and error range using ANN, allowing the controller to 
successfully track the desired trajectory. The amelioration was noticeable on the position tracking where the 
controller became insensitive, and the disturbance had no effect on position tracking. 


— Estimated value of M 0.5 — Estimated value of D | 
6f -- Real value of M ] -- Real value of D 
5 _ 0.4 
oH 
hd 
m4 ig l P 
oP a = 0.3 
= 2 
panty, e 
2 4 A 
1} 0.1 
0 i : 5 0 
0 1 2 3 4 5 6 7 8 0 l 2 3 4 5 6 7 8 
Time (s) Time (s) 
(a) (b) 
i 
FS RE PE PE mn vo, | 


Y 


= - Estimated L without ANN 
= Estimated L with ANN 
*. Real value of L 


(c) 
Figure 6. Parameters estimates: (a) mover mass, (b) friction coefficient, and (c) load force disturbance 


4. CONCLUSION 

This paper has proposed an adaptive backstepping control structure using artificial neural network 
for load disturbance estimation, to control a LIM position. The controller’s structure was designed and 
parameter update laws were extracted, to track a desired position trajectory, assuming that some LIM 
parameters are uknown. The robustness of the proposed controller was checked under different operating 
conditions. The obtained simulation results have shown a great position tracking performance compared to 
conventional backstepping. And they have proven it’s insensitivity against parameter variations, and it’s 
ability to compensate load disturbance effect without overshoot, thanks to parameter update laws and ANN 
estimation structure. Finally, the research directions of this paper can be extended in the future, to the use of 
ANN backstepping controller in tracking the speed of LIM. Which means considering the end effect 
phenomenon when modeling LIM. 
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